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^ ■ Abstract 



We discuss a relativistic theory of the atomic nuclei in the frame- 
work of the hamiltonian formalism and of the mesonic model of the 
nucleus. Attention is paid to the translational invariance of the theory. 
Our approach is centered on the concept of spectral amplitude, a func- 
tion in the Dirac spinor space. We derive a Lorentz covariant equa- 
tion for the latter, which requires as an input the baryon self-energy. 
For this we either postulate the most general Lorentz-Poincare invari- 
ant expression or perform a calculation via a Bethe-Salpeter equation 
starting from a nucleon- nucleus interaction. We discuss the features of 
the nuclear spectrum obtained in the first instance. Finally the gen- 
eral constraints the self-energy should satisfy because of analyticity 
and Poincare covariance are discussed. 
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1 Introduction 



The purpose of this paper, which extends and deepens a previous one[|T|, is to 
outhne a theory of finite nuclei cast in a Lorentz and translationally invariant 
form. 

Extensive work has been carried out in the past to describe a nucleus 
in a relativistic regime (see for instance 0, |^ for comprehensive reviews). 
This has been done essentially in the framework of the Dirac's theory, either 
phenomenologicallylQ, |, 0, i.e. by suitably parametrizing the self-energy felt 
by a relativistic nucleon when interacting with a nucleus, or microscopically, 
i.e. by evaluating the self-energy starting from a mesonic theory, like the 
Bonn potential g). 

However, the Dirac's theory, as it is usually applied, breaks translational 
invariance. This appears to be serious at high energies as already indicated 
by nonrelativistic calculations, which point to the growing importance of the 
center-of-mass component in the nuclear wave function when the momentum 
transfer becomes large0. 

A notable extension of the Dirac's theory is represented by the Quantum 
Hadrodynamics (QHD) of Waleckafl^, |lT|, a quantum field theory model 
based on a lagrangian which describes nucleons interacting via the exchange 
of isoscalar scalar (a) and vector {uj) meson fields. 

Notably, QHD has been remarkably successful in dealing both with nu- 
clear matter and finite nuclei in the so-called mean field approximation 
(MFA), where one replaces the meson field operators by their ground state 
expectation values ignoring furthermore the negative energy states. However 
when one goes beyond a MFA to include correlations among nucleons (via, 
e.g., a Bethe-Salpeter-type equation or truly quantum field theory effects, 
as those stemming from vacuum polarization) then QHD faces severe prob- 
lems. Indeed, although great care is paid in its framework to the problem 
of renormalization, the loop expansion appears to be, to say the least, quite 
poorly convergent. For example the two loop contribution to the binding en- 
ergy turns out to be very large at the nuclear matter density[|l^] because of 
the unsatisfactory ultraviolet behaviour of the theory, a well known feature 
common to all the non asymptotically free model. One can of course deal 
with these difficulties by introducing ad hoc form factors to cut off the large 
momenta entering the diagrams needed to express the physical observables of 
concern, but in so doing one spoils QHD of its character of a truly hadronic 
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field theory. 

Our approach differs from QHD in several respects. First it is centered on 
the concept of single particle orbital. In fact, notwithstanding the difficulty 
of defining a wave function, in particular a single particle one, in a relativistic 
framework, yet it is still advantageous to deal with the nucleus in terms of 
a relativistic extension of the quasiparticles, the building blocks of atomic 
nuclei, as introduced by Landau and then by Migdal[14|. 

To achieve this goal it is crucial to introduce the so called spectral ampli- 
tude (pi, to be later defined, a suitable tool for expressing the nuclear binding 
energy, the nucleon momentum distribution and therefore the cross sections 
for a variety of nuclear reactions (for example the exclusive {e,e'p) one ||15||) 
in a Poincare invariant form. 

A second difference from QHD arises from the recognition that the con- 
nections between 0; and the observables above referred to are more naturally 
established in a hamiltonian framework where they can be cast in the form 
of exact relations. Of course the problem of the infinities arising from the 
presence of the Dirac's sea exists and we deal with it empirically by means 
of ad hoc subtractions. 

Moreover the hamiltonian scheme has the advantage of allowing us to 
fulfill more easily (at least formally) Poincare invariance when treating finite 
nuclei. As a notable consequence the Dyson's equation for the fermion prop- 
agator, from which an equation for 0/ is derived, exhibits an algebraic, rather 
than an integral, structure even for finite nuclei. Obviously here the difficulty 
is met of properly defining and calculating the self-energy of a baryon or of 
a meson (see sect. §). 

In this connection we adopt a pragmatic view, namely we write the most 
general Poincare invariant self-energy in terms of functions meant to be fixed 
by the phenomenology 0, |16|, |, and we explore the consequences it entails 



on the spectrum of the system. 

We also investigate the alternative option of calculating the self-energy: 
here our approach is based on a Bethe-Salpeter-type equation^], with a 
phenomenological nucleon- nucleus interaction as an input. We thus account, 
respecting covariance, for the repeated interactions a particle undergoes in 
its way-out from the nucleus (final state interaction). It is worth noticing 
that this approach has the merit of transparently displaying for a convenient 
choice of the interaction and for an infinite homogeneous system how our 
hamiltonian theory reduces to QHD at the mean field level. 
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This paper is organized as follows: in sect. ^ we outline our hamiltonian 
framework, introduce the spectral amplitude and establish the connection of 
the latter with various physical observables. In sect. ^ we discuss a simple 
example to illustrate how mesonic degrees of freedom can be accounted for 
by means of a canonical transformation. In sect. |^ we study the algebraic 
Poincare-covariant Dyson's equation for the fermion propagator in a finite 
system. The Lehmann representation for the propagator is also analyzed. We 
study then the general structure of the self-energy, its link with the final state 
interaction via a Bethe-Salpeter equation and its analytic properties. Finally 
we shortly analyze the resulting energy spectrum of the nuclear system. 



2 The formalism 

Let us consider a system of nucleons interacting via the exchange of mesons. 
For simplicity we confine ourselves in this section to consider the exchange 
of a scalar meson. Euristically the latter simulates the exchange of a pair 
of pions, a microscopic process viewed as the main source of the nuclear 
binding. Clearly our model should then be generalized to more realistic 
situations encompassing several different mesonic fields to account, e.g., for 
the short range repulsion among nucleons. 
Its hamiltonian reads 



where 



H = H% + Hl + H' (1) 

Hi = Jd'y^iy){-t-fV + m)^P{y) (2) 
= \ j d'y [a\y) + (Vo-(y))2 + mW{y)\ (3) 
H' = gJd'yi^{y)aiy)^Piy), (4) 

ip and a being the fermion and meson field. 

A flaw of (|l]) stems from the absence of a self- interacting term Xa'^, needed 
to bound from below the spectrum of (|ip. However since, as we shall see in 
the following, the stability of the a field is here ensured by the presence 
of the nuclear medium, the renormalized value of A can be assumed to be 
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positive and small: hence the Acr^ term is not expected to significantly affect 
the predictions of the theory and can thus be safely omitted. 

A further shortcoming of (|I|) is connected with the difficulty of defining 
a vacuum for a bound, finite, relativistic Fermi system. Indeed the standard 
expansion for a fermion field, namely 



^(y) = E 



m r 



(27r)3 ek 



iky 



(5) 



{a^s and bks being the free nucleon and antinucleon annihilation operators 
and ek = Vk^ + m"^), is not appropriate, being incompatible with the spatial 
confinement of the atomic nucleus. 

However, since we do not pretend to obtain here a solution of the problem 
of confined states in QFT (in fact we rather search for relations between 
physical observables and the spectral amplitude) we shall still quantize the 
fermion field according to (^, but without introducing normal products at 
all. We shall instead resort to the introduction of compensating terms, in the 
form of vacuum subtractions, in order to avoid the occurrence of infinities. 

On the other hand, as far as the mesons are concerned, we assume their 
coupling with the fermion field to be so weak to prevent the occurrence of 
bound meson states. Therefore the expansion 

d^k 1 



Cke'^-y + cle 



iky 



(6) 



(cfc and c^k being the free meson annihilation and creation operators and 

and 



y k^ + m^) is likely to be still warranted. 
In line with the above considerations we shall rewrite equations 
as follows: 



d^k 



H' 



2u;fc(27r 
dH: 
2uJk{2n 



;^ ^kC'kCk 



3 Apk + pW 



where 



Pk = g d^yilj{y)ip{y)e 



-iky 



(7) 
(8) 

(9) 
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Our choice of the normahzation of the fields is such that the commutator 



[cfc,4,] = (27r)32^fe5(k-k') , (10) 
holds for the boson and the anticommutator 

{aL,afcv} = (27r)3^(5,,,(5(k-k') (11) 

for the fermion fields. 

Let us now introduce the nuclear states |p, {A)m >, identified by the total 
3-momentum p, baryonic number A and by a set of internal quantum num- 
bers denoted by m. With the notation |p, A > we refer to the internal ground 
state. Since the relativistic normalization of a state is statistic-dependent, 
we have to choose the parity of A. Let us assume A to be even. Then the 
normalization reads 

< p, A\p', A >= {27rf2po6{p - p') . (12) 
Of course for the ground state 

P^\p,A>=p^\p,A> (13) 

and 



< P,A\H\ 


\p,A> 


< P,A\ 


\p,A> 



Po= ^' ' , ' , <0\H\0>= JMl + p^ . (14) 



Clearly < p,yl|if|p, A > diverges owing to the contribution to the energy 
arising from the Dirac sea, hence the vacuum subtraction in ([T^). 

In parallel to (|T3D for a (A — 1) system, left in a quantum state identified 
by the index /, we have 

P,\p,iA-l)i>=p'^\p,iA-l)i> , (15) 

where 

pi = (pI p) (16) 
and 



p[,^6f-i(p) = V(M^i)2 + p2, (17) 
but now the fermionic normalization 

< p, (A - l),|p', (A - 1), >= (2nf-tl^Su>5ip - p') (18) 
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holds. 

Let us now introduce in the Dirac spinor space the function 



(l^ikpiy) = Ciik, p) < p - k, (A - |p, A > (19) 

with 



^ 2poer(p-k) ' 

the field ipiv) being now in Heisenberg picture. The above is the key quantity 
of our formalism, in terms of which we shall express the ground state energy 
and, in a parallel research, the structure functions of the nucleus. In sect. ^ 
we shall derive the equation obeyed by (pi. 

The amplitude (pikpiy) is the tool we need to generalize the quasipar- 
ticle concept to a relativistic regime: in fact it can be shown that in the 
nonrelativistic limit and in the single-particle approximation (pi exactly re- 
duces to a one particle wave function. Furthermore the coefficient Ci(k, p) is 
such to account for the difference between the relativistic and nonrelativistic 
normalization of the states. 

Of course the spatial dependence of (pikp{y) is fixed by translation invari- 
ance. Indeed, from 

^(y) = e'^-y^{Q)e-'^-y (21) 



one gets 
where 
and 



(t>ikp{y) = m,v)e-'''-' (22) 
0z(k, p) = a(k, p) < p - k, (A - l)dV'(0)|p, A > (23) 

^ (e«(k,p),k) (24) 
ei{k,p) = p,-ef-\p-k). (25) 

Owing to the definition ([T9|) 0i(k, p) transforms under the Lorentz group as 
a Dirac spinor. Moreover, it can only depend upon the Lorentz-invariant 
quantities p^, (fc')^ and p ■ k\ but the mass shell conditions (p^ ) and (p5|) 
leave only the 3-vectors p and k as independent variables. 
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The quantity |0i(k, p)p is the probabihty density (homogeneous in space 
and time) of destroying a nucleon (or creating an antinucleon) with 4- 
momentum ki in the nuclear ground state leaving the residual nucleus in 
the quantum state |p — k, (A — 1); >. Notice also that eqs. (^) to (pSj) hold 
valid in the nonrelativistic limit, where the center-of-mass can be correctly 
defined and its motion accounted for. If however the center-of-mass is not 
properly singled out, as is the case, e.g., of the nuclear shell model, then 
obviously the invar iance for space translations is lost. 

Let us now illustrate how the energy of the nucleus and the momentum 
distribution of its constituents can be expressed in terms of 0/(k, p). 

To start with, notice that 

i:/|^i*(k.P)r- p6) 



J d'y <p,A\^\y)^{y)\p,A>=A + 4Q J 



where the term proportional to the normalization volume Q arises from the 
Dirac sea and would be absent had we employed a suitably defined normal 
product. 

Furthermore, from (EH) and using closure, one gets: 



E/'(0j^Kk,p)|0z(k,p)r (27) 



< P,^|P,^ > 



J d'y <p,A\i;\y)ipo-H)ijiy)\p,A> 



which, exploiting the commutation relations for the field ip, can be recast in 
the form 

E / ^eKk,p)|0Kk,p)r (28) 

= 7 T^^<P,^|(i^^ + i^')lP'^> ■ 

<p,A\p,A> V y 

The above relation expresses the energy of the system associated with the 
fermionic degrees of freedom. This quantity can of course be experimen- 
tally inferred from the nuclear separation energy. The "total" energy of the 
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nucleus, which is also experimentally known, reads instead 

Po = E / 7^^/(k,p)|0Kk,p)|' (29) 



^ r <fk <Y>,A\ukclck\v,A> ^ ^ 
J 2u}k{2TiY < p, A|p, A > 

One thus sees that in principle the possibility is offered, within our hamilto- 
nian model, to estimate the contribution to po arising from the mesonic com- 
ponents in the nuclear ground state. This is expected to be small: should not 
this be the case then the mesonic field expansion (|^) would not be compatible 
with the normal product introduced in (|^). 

Let us now somewhat generalize equation p9|). For this purpose we con- 
sider again the baryon contribution to the total energy of the nucleus, intro- 
ducing the analogous of the spectral function in the nonrelativistic theory, 
namely 

S{k, p) = 27r E |0,(k, p) |2 5(A;o - e^k, p)) - S,^,{k) (30) 

which, as it is well known, is connected to the correlation function via a 
Fourier transform 

S{k,Y>) = T^^Tr / dS < p, A|V^t(^)^(o)|p, A > e-'^-y - S^,,{k) (31) 



and is crucial for expressing the (e, e'p) phenomenology [jl5 



In term of the above, and reminding the need of subtracting the vacuum 
contribution 

S,,,{k) = %-Kn5{ko + eu) , (32) 
the baryonic part of the energy can then be expressed as follows: 



nucl ' " 



pT = J (^^o5(A:, p) . (33) 
The mesons contribution can be similarly handled by setting 

S'Ak, P) = 2vr^Pl^M^L^5(A:„ _ (34) 
2uJk < p,A\p,A > 
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in terms of which the mesonic energy reads 



Pr = /(0^o5°(A:,p) . (35) 

The previous result (^) may then be generahzed by considering the total 
4-momentum, whose expression 



(27r) 



immediately follows from the above definitions. In ( P5| ) S and can clearly 
be interpreted as the 4-momentum distributions for nucleons and mesons 
respectively. 

Note the asymmetric nature of (p^): S{k,p), as defined in eq. (pO|), 
embodies the exact spectrum of the system, whereas S^{k,p) contains the 
free meson spectrum. Thus S{k, p) provides all the needed informations 
about the interacting system: what is left out is the hamiltonian of the 
free mesons. Of course one could have started by defining an exact spectral 
function for the mesons, which would then contain the full spectrum of the 
system: then a free spectral function for the nucleons would be left out. 

An energy integration leads to the 3-momentum distribution of the 
fermions 

n(k,p)=5:i0,(k,p)r (37) 

I 

and of the mesons 

2uk < p, A\p, A > 
respectively, so that the total momentum p of the nucleus will be 

/d'^ k 
^kMk,p)+n.(k,p)} . (39) 

Finally, note that S, and S'vac are Lorentz-covariant and the normal- 
ization of S is 
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Clearly 

provides the average number of mesons inside the nucleus. 



3 A Simple Case 

In the previous section we have outlined a Lorentz-invariant scheme for 
treating the atomic nucleus. In practical cases however one is often 
forced to resort to a single-particle approximation (Dirac-Hartree-Fock 
approximation [0, 0). In so doing translational invariance is generally lost 
and the mesonic degrees of freedom are dropped out of the problem. 

Here we examine, in a simple case, how a description of the nucleus solely 
in terms of nucleonic degrees of freedom can be worked out respecting Lorentz 
and translational invariance. 



For this purpose we start by writing the purely nucleonic state vector [IS 
(sometimes also called Fock's line) 



\p,A>F={2nr^^ (42) 

n (2^~^P^^^' ■ ■ ■ ' ^^^^^^ • • • «LlO > ^(P - E ki) 

where k; = (kj, Sj), Si embodies the spin-isospin variables, a^. and a^^ are the 
nucleon creation and annihilation operators, ek = \/k^\-rn? and the vacuum 
definition 

aks\0 >= bks\0 >= (43) 

is assumed to hold. 

With the normalization (0) for the states \p, A >, and owing to (p^), 
the normalization of reads 



^'^^ '^\Mh,...,kA)f6{p-j:k,) = l (44) 



A(2vr)^6,, 

where ipp(k) is also fully antisymmetric in its arguments. 
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The advantage of ( ^21) is of rendering straightforward the connection with 
the non relativistic theory: indeed to recover the latter it is sufficient to 
replace ipp{ki, . . . , fc^) with, e.g., a shell model wave function in momentum 
space. 

In addition (0) allows us to express the energy of the nucleus and the 
4-momentum distribution of the nucleons in terms of the 3-momentum dis- 
tribution only. In fact since the 3-momentum distribution reads 

^Fn. _^ F< p,A\al^aks\p,A>F 

" ,<p,A\p,A>, ^^^^ 

{s being the spin index), the energy will be given by 

p _ F<P,A\J (0gefcaLafcs|p, A >f 
^° ~ F <P,A\p,A>F 

= E/(0^.<(k,p). (46) 



Note that the Lorentz-invariant 3-momentum distribution (^51) is normalized 
according to 

as it easily follows from (0), and that the 4-momentum distribution is ex- 
pressed in terms of the latter as follows 

S^ik, p) = 2ttY, nf (k, p)5iko - e^) + Sn^diko + e^) (48) 

s 

(of course 5° = 0). 

Also, for sake of completeness, we write down the expressions for and 
in terms of the function ipp{k): 



n 7:711- \Mks, k,,..., kA)\' 6{p-k-Y: k,) (49) 

i=2 \'^'^) j=2 

= (2^)' / n tC^- mk,, ...M)? e,AP - E k,) . (50) 
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The problem with ( ^21) is that it can never correspond to a bound state. 
Indeed, in the rest frame, the energy of a bound system clearly satisfies the 
condition 

e < mA . (51) 
On the other hand, from (^) and (^7|) it follows 



Po 



7^^^<(k' P) > E / 7^<(k' p)=mA. (52) 



(2vr)3 - (2 



The above result is obvious because the mesons do not intervene when the 
expectation value of the hamiltonian is taken in the state (|4^) . However it 
points to the need of inserting the mesons in the model, preserving, possibly, 
a structure like (^) for the nuclear state. To illustrate how this scope can 
be achieved in a simple scheme, let us consider an infinite Fermi gas at zero 
temperature in the rest frame. 

Notice that the well known momentum distribution for such a system, 
namely 

(k) = ne{kF - k) , (53) 

is not Lorentz-covariant. Yet we are willing to sacrifice temporarily Lorentz 
covariance in order to take advantage of the transparency of the Fermi gas 
model. 

We start then from the equations of motion |TP|, |2D[: 

(□-"^a)^ci = g<i^i^> (54) 

(i7 ■ (9 - m - gcTci) ^ = . (55) 
These, for a uniform system (the Fermi gas), simplify to 

(Tci = — ^ Kipip > , (56) 

{i-fd-m)ip = , (57) 

9^ - 

m = m < ipip > . (58) 

Thus an effective mass, always lower than the bare one no matter which sign 
the coupling constant has, naturally appears. 



13 



We then perform the canonical transformation 



a ^ a — (Tc\ (59) 

to get rid of the average value cTci: this amounts to replace the free nucleons 
and antinucleons with dressed ones. These may be viewed as baryons carrying 
a cloud of (T-mesons, but otherwise unaffected by the interaction according 
to eq. (p?!). To evaluate the energy of the system we accordingly use the 
field representation 

^(y) = E / yf^- (n..a,.e^'=-^ + hsPle-^''-^) (60) 

where the operators a and /5 create and annihilate dressed nucleons and 
antinucleons (e^ = A/k^ + m^) . Proceeding now along the previously outlined 
scheme, but utilizing the representation (|60|), which entails a vacuum different 
from (|3) and, consequently, a different subtraction, we get 



1 ^ f d^k ^ 



, J (Znr ek 



<p,A\p,A>^J (27r)3 gfc 
+ is^mM-4n/|^jf,. (61) 

Here the first term has the lower bound rfiA < mA, whereas the second one 
arises from the classical meson field via the canonical transformation (|5P|). 



Thus a bound nuclear state becomes possible, depending upon the value of 
the coupling constant g and of the nuclear density. The last term in (^) is 
again due to the lack of a suitably defined normal product and its structure 
reflects the presence of two vacua, one associated with the free particles and 
the other with the dressed ones. 



4 General case 

In this section we shall examine how a many-body theory for a finite system 
can be cast in a relativistic framework still preserving Poincare invariance. 
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4.1 The Single-particle Green's Function 

We start from the single-particle Green's function 

. <V,A\T{^{y),^\y')]\Y^,A> 
y ^ = <p,A|p,A> (^2) 

which obeys the Dyson's equation 

{zj-d- m) Gp{y, y') + j d^z^^^y, z)G^{z, y') = ^^5\y - y') . (63) 

Since we employ momentum eigenstates Gp can only depend upon y' ~y 
(translational invariance). It will carry however an index p representing the 
total momentum of the nucleus (the same applies to the self-energy Sp as 
well). Accordingly, we can work in the 4- momentum space, where eq. (|6^) 
becomes 

i-fk-m-Ep{k))Gpik)=jo . (64) 
Obviously, eqs. (|63D and (|6^) are meaningful only if a prescription is given 



for calculating the self-energy. This is far from being trivial in a relativistic 
context. A perturbative approach will not work, since a bound state of a 
finite assembly of nucleons cannot be obtained from a non interacting system 
via the summation of Feynman diagrams. Rather a collection of nucleons 
bound in an external (possibly self-consistent) well (as in refs. [T7[| ) 
might be taken as a starting point. This however entails to give up Lorentz 
and translational invariance from the very beginning in the hope that the 
summation of appropriate classes of perturbative diagrams will restore those 
symmetries. If this turns out to be the case, then eqs. (^) and (|6^ are 
meaningful. 

Here our standpoint is that Sp exists, is calculable, Lorentz-covariant and 
transforms like a 7-matrix. As a consequence Gp will exist as well. However 
to write down a Dyson's equation of the type 

Gp = Gp + GpTipGp , (65) 

preserving at the same time the Poincare invariance together with the exis- 
tence of a bound nuclear state, turns out to be a task beyond reach. 
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Therefore, we start from the spectral representation of Gp itself, namely 



< p, A|p, A > fco - £^/(k, p) - ir] I h- e'Kk, p) + ^^7 j 

(66) 

which introduces, in parallel to the previously discussed 0i(k, p), a new func- 
tion 

X/(k, p) = G(k, p) < p + k, (A + l)d^t(0)|p, A > , (67) 

with 



Q(k,p) 



(68) 



whose modulus square yields the probability of creating a nucleon (or de- 
stroying an antinucleon) with 4-momentum ki on the nuclear ground state 
|p, A > leaving the residual system in the quantum state |p + k, (A + 1); >. 
Moreover, ef'*'^(p + k) and e;(k, p) refer to an excited state with baryonic 
number A + 1 and are defined in analogy with the A — 1 case. 

It helps, in grasping the significance of the singularities in the propagator, 
to think to the infinite Fermi gas where one has 

G^^oik) = + ^ [ ^(^ - (69) 
^~ 2v^P~+^r? [ko — \^k^~^^rn? + ir] 

^ e{kF-k) 1 ] 



/cq — yl^^rrri? — ir] k^ + yli^^+rri? — it] j 

Here the antinucleon states are entirely embodied in the last term of the 
rhs. This separation however becomes impossible when the interaction is 
introduced. Accordingly in ( p6| ) the antinucleons are actually embodied in 
both terms of the rhs (not only in the first). 

To better understand the above spectral representation it might help to 
think to a system in which either a nucleon has been annihilated or an antin- 
ucleon has been created. This picture strictly speaking is not correct, since 
the nucleus is actually described by a superposition of states with different 
number of nucleons and antinucleons; we can however start by considering an 
assembly of noninteracting nucleons in the rest frame and then annihilate a 
nucleon or create an antinucleon. In this way two sets of states are obtained, 
denoted, with self-explaining notation, {A — IN) and {A + IN) respectively. 



16 



Then switching on the interaction we generate two classes of physical states, 
still somewhat linked to the original {A — IN) and {A + IN) ones. The pos- 
itive energy solutions are associated to the former, the negative energy ones 
to the latter, both stemming from the singularities of the first term on the 
rhs of (^). The singularities of the second term relate to {A + lA^) states, 
i.e. states in which a A^A^ pair exists. 
Of course 

^+4/^. (70) 



n J {2tt) 

By inserting then ( ]66| ) in (0) one obtains the following eigenvalue equa- 
tion for 0i(k, p) in terms of the a, (3 representation for the Dirac's matrices: 

Gpi(^)0Kk,p)= (71) 
= fe^k, p) - a . k - /3m - /3Sp(fc) ) 0;(k, p) = . 

An analogous equation holds as well for the x^k, p)- 

In (|7TD the eigenvalue £i(k, p), as told by eq. (|25|), represents the dif- 
ference between the ground state energy of the system in the state |p,A > 
(moving with total momentum p) and the energy of the system in the state 
|p — k, (A — 1); > (moving with momentum p — k and intrinsically left in 
the quantum state /). It can be regarded as the energy of the quasi- hole. On 
the same footing e;(k, p) represents the energy of a quasi-particle. 



4.2 The Self-Energy and the Final State Interaction 

The crucial quantity in the many-body equation for the spectral amplitude 
is clearly the self-energy. In this subsection we shall obtain the latter from 
the interaction between a nucleon and the daughter nucleus, thus linking the 
self-energy to the so-called final-state interaction. An attempt in the same 
direction has been performed in ref. 0]. 

Let us introduce a new field operator '^{Y) associated with the A — 1 
nucleus viewed as a composite fermion which can exist in different internal 
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states: 

(72) 

Here the operator Aip annihilates a nucleus of baryonic number A — 1 in the 
excited state / and with total 3-momentum p (we neglect, as irrelevant in the 
energy range of actual interest, the component of the field that creates the 
corresponding antinucleus). Uip{Y) is a Dirac spinor, normalized in Fourier 
transforms according to 

UipiY) = [/^pe^^P 

UipU.p = '-^^Su> . (73) 

Next we introduce the four-point Green's function 

_ <P,A\T{iiymY»Ky'mY')}\p,A> 
y^[;ll,r,y,r )- <p,^|p,A> ^ ' 

which obeys the following Bethe-Salpeter-like equation 

gp{y, Y; y\ Y') = gl{y, Y- y\ Y') - i j d'yi d% d'Y, d^Y^ 
xgl{y, Y- 2/1, ri)Vp(2/i, Fi; 2/2, ^2)6^^(2/2, Y<,- y\ Y') (75) 

where is the propagator of the A^^ nucleon and of the {A — 1) daughter 
nucleus in the absence of their mutual interaction Vp. 

Since all the functions Qp, and Vp carry a definite total momentum p, 
they are best handled in Fourier transforms: accordingly let us set 

J^{k,k';p) = J d^y dy d'^Y d^Y' (76) 

where T can be Q ^ or V. In this representation reads 

Q\k,l^-,V) = {2'K)%k-k')g\k)G\_,{p-k) . (77) 

Here g^{k) is the relativistic free nucleon propagator in the vacuum (not to 
be confused, of course, with the G° appearing in eq. (|65|) ). Indeed the latter 
is ill-defined, whereas the former (disregarding renormalization effects) reads 

= fco - a ■ k - /3m . (78) 



18 



The relativistic free propagator for the {A — 1) nucleus (not interacting with 
the A^^ nucleon), is given by 

I (p)Po-e/ (p) + 
The Fourier transformed Bethe-Salpeter equation reads then 

gik,k';p) = {27:fg\k)G\_,{p-k)5{k-k') (80) 

+g'^{k)G\_^{p-k) I d^kV{k,k;p)g{k, k'-p) 

We introduce now the Bethe-Salpeter wave function in line with the usual 
convention ^ ^| 



(^p(y,F) = < 0|T{V^(y)vI/(F)} |p, A >= j ^^p(fc)e-'=^-(^-^)^ (81) 

whose spectral representation reads 
^p(fc) = (82) 

ko - H + trj ko - po + H - tr] 



P being the 3-momentum operator and Po = y M\ + both in eqs. (^) 
and (P^). 

From (|80D an equation for v^p(/i;) is then deduced by writing down the 
Lehmann representation of g{k, k';p) and selecting those contributions hav- 
ing the vacuum as intermediate state. This scope is achieved by integrating 



over po around a small circle surrounding the point Po = yp"^ + M\ which 
indeed corresponds to the pole of the propagator having the vacuum as an 
intermediate state. Thus we find 

j^g{k,k'-p) = ^l{k)v^{k') (83) 

and since is regular at po = y^p^ + M\ one gets, by integrating eq. (^) 
over pq, 

Vp{k) = -zf^/{k)GUiP-k) J 0^^Vik,k';p)^^{k') . (84) 
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Remarkably, the Bethe-Salpeter wave function (pp{k) turns out to be con- 
nected with the spectral amplitude 0i(k, p). In fact, as seen from eq. (|8^) , 
(Pp{k) is a meromorphic function in the variable fco, with poles located at the 
energies of the free baryons (first term of the rhs of eq. (0) ) and at the 
energies of the quasiparticle (second term). We can thus insert a complete 
set of eigenstates of the {A — 1) system in the second term of the rhs of eq. 
(82), integrate over a small circle around the selected pole and, by reminding 
the definition (ESl) of we find 



From the above, using (^) and (|82D, the connection we are searching for, 
namely 

MK P) = Q{k, p) / ^<_,)V^p(A;) (86) 

follows. 

Next we act on (|8^) with [(7"]"^ and again select the poles by integrating 
over fco around the l^^ eigenvalue of the [A — 1) system getting 

(si (k, p) - a ■ k - (3m) (k, p) = (87) 

a(k,p) / 0^^u{^_^^V{k,k';p)vp{k') . 

The above equation provides a tool for calculating microscopically the 
self-energy: indeed given an effective nucleon-nucleus interaction V{k, k'; p) 
the Bethe-Salpeter wave function (pp{k) can be found and the product of the 
baryon self-energy with the spectral amplitude (/);(k, p) can be determined. 

Note also that the finite size of the system is actually embedded in 
V{k,k';p). 

Clearly for V{k,k';p) a model is needed. Interestingly, if we choose the 
expression 

V{k, k'- p) = i2nfS{k - k') (pa + cu) (88) 

for the effective interaction, then the well known self-energy of the QHD, 
namely 

Sp(fc) = a + f3uj , (89) 
is recovered, a and u being the standard scalar and vector fields of QHD. 
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4.3 The Structure of the Self-Energy 

Having examined a microscopical model for the self-energy and its link with 
the final state interaction, we now write down its general expression, namely 

E(fc2, k-p) = S + ^,V^ + ^ [7;., 7.] T^" , (90) 

which fulfills the constraints imposed by Poincare covariance. Our purpose is 
to explore the predictions that can be made on the relativistic quasiparticle 
spectrum without specifying the details of the baryon self-energy. In (pO|) 



V = ki'V + p^'V (91) 

and 

T^^ = k^'p'^T (92) 

are a Lorentz vector and a second rank tensor respectively, whereas S, V and 
V are scalars, in general complex, whose dependence is upon the invariants 
k'^ and k ■ p. Alternatively, we can choose as independent variables, fc^ and 
{p — kY". this might be more convenient as it will be discussed in sect. ^3 . 
Our expression for Sp is similar to the one of ref . |]16| , often considered in 



the literature, but we remind that translational invariance strongly constrains 
its structure. We also remind that terms with 75 are neglected as parity- 
violating; concerning the impact of other discrete symmetries on Sp we refer 



the reader to [16 



With the self-energy (^OD the inverse propagator entering eq. ([7l|) reads 
G~^\k) = {l-V) - t;po - « • (k - - (3m* - [/k, /p]} (93) 



where the definitions 



and 



v 

V = j^, (95) 
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have been introduced. The above Green's function Gp describes a moving 
quasiparticle interacting with the daughter nucleus. However, before dis- 
cussing the physical significance of the functions introduced so far, namely 
m*, t and V , we first need to find the poles of the Green's function, i.e. 
to solve the eigenvalue equation for the spectral amplitude 0^, and then to 
study the analytical properties of Sp. 

4.4 Solution of the Eigenvalue equation 

To solve eq. (^) we look for the poles of the propagator 

^^^^^ " (1 - V) (P - mf + 17^) ^^^^ 
jfco -vpo-(x■{\^- VY)) + /3m* - [/fc, /p] /jj 

where a new mass 

m\ = m*^ + 2vk-p- t; V + | A;^ - ( '-^-^ ] I (gg) 




has been introduced. The small positive imaginary part in the denominator 
is meant to be effective as far as m*, v and t are real. 

Clearly the poles of Gp coincides with the eigenvalues e/(k, p) of (|7l|) and 
are the roots of the equation 

P = mjik'^ , k ■ p) . (99) 

We shall discuss the regularity of the functions m*, v and t later on. In the 
regions where they are well behaved and m\ is positive definite the nonlinear 
equation (B^) might display positive 



+ ^lk=ei (100) 

and negative 

er = -^Jk + ml 1,,,^^- (101) 
solutions. Notice that these in the rest frame become 



= MAVi±J{l-tf)k^ + mf (102) 
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where the suffix / reminds that the various quantities have to be evaluated 
at ko = ef. 

We thus see that vi yields an overall frequency shift, whereas m*i sets 
the splitting between the two branches of the spectrum. Noteworthy is that 
an energy shift, a mass renormalization and the introduction of an effective 
momentum 

keff = Vl -t^k (103) 

allow a recovering of the structure of the Dirac's equation eigenvalues. In 
particular the scalar component of Ep renormalizes the mass, the vector 
one partly renormalizes the Green's function residua and partly shifts the 
spectrum and finally the tensor term modifies the momentum. 

Of course m\ is not an even function of fco (due to its dependence upon 
p ■ k), so that in general e'^ ^ ■ Also, the existence of one solution does 
not imply, in general, the existence of the other one. 

From the analogy with the Dirac's equation we can interpret these so- 
lutions as describing the propagation of a nucleon (or of a hole, depending 
upon the boundary conditions) or of an antinucleon in the medium. 

Finally we note that the positive branch is the meaningful one in the non 
relativistic limit, where 

ei{\i)^ m* + MaVi + '^^'}^^ ■ (104) 

2m* 

Now we search for the spinor structure of the spectral amplitude 0/. The 
equation (^) can be recast, with the help of (|97|), in the form 

(^£i{k,p) - vipo + cx{k - vip) + (3m* - [/fc |fcg=e,(k,p), /p] /?} 0Kk, p) = 

(105) 

which defines </)i(k, p) up to a normalization. 
Now for positive frequency solutions we get 

^r(k,p)=i:(i + ^*)"^j(k.P)L„^,. (106) 

where 

".-(/^f.S.). (107) 
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K = k(i- e 



M = m* - 



( m* k ■ p 

P \V TTTrt 



M 



A 



Ml 



Ma 



k ■ p — vp" 



and 



For negative frequency solutions we have instead 



ko=£ 



with 



IS + M 
2M 



(108) 
(109) 

(110) 
(111) 

(112) 



In (|107|) and (112) <^a is the Pauh spinor corresponding to the spin polariza- 
tion s. We normalize the spinors u and v according to mm = 1, vv = — 1, 
while the normalization coefficient A/" remains to be fixed. This is done by 
observing that near the pole £;^^(k, p) the Green's function behaves like 



G ~ 



but on the other hand 



2po 



^>f(k,p)0±^(k,p 



< pA|pA > ko — £:(k, p) — irj 



Gr 



{ko-ct-k-pm- 



-1 



By comparing ( [113|) and ( |114| ) near the pole the relation 



2po 



< pA\pA > 



(k,p) 



VPO 
kp 



(113) 



(114) 



(115) 



ko=£. 



which fixes the normalization of 



'I 5 



follows. 



4.5 Properties of the Self-Energy 

In this subsection we first shortly study the analytic properties of the self- 
energy starting from its microscopic definition. For this purpose we consider 
the equation of motion for the field operator ip{y) 

{i-i ■ d -m-U)i){y) =Q (116) 
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where the operator lA is defined through the equal time commutator 



W){y) = [i;{y),H\yo)] . (117) 

The operator U here introduced is model dependent. For our hamiltonian 
(|T]) it reads 

iU^|J){y) = ga{y)i;{y) . (118) 

Now, by operating with the operator on the Ihs of ( |116| ) on the Green's 
function definition and by exploiting the Dyson's equation (|63D, one obtains 



(SG) (y, y') ^ I rfSS(y, z)Giz, y') (119) 

<Ap|T{(W^) {y),l^\y')]\A^> 

= ^Wo ^ T, : • 

<p,A\p,A> 

From ( [L19| ) the Lehmann representation 

< A\A > ^ ' 

< A|^t(o)| - k, (A - 1); X -k, {A - l)i\{Uijm\A > 
I ko - eiik) - ir] 

^ ^ <A\{utpm\- k, {A + i)i X -k, {A - iwmA > \ 




ko - e(k) + if] 



is deduced in the rest frame. 

Let us now restrict ourselves to the discrete spectrum only. Here, as a 
function of ko, G{k) (which, strictly speaking, is a distribution) displays a 
pole in correspondence of each discrete eigenvalue. Owing to (|120|) so does 
SG, hence E must be regular. Furthermore a zero of G certainly exists be- 
tween two neighboring poles, where instead EG is regular. Hence S must 
diverge where G vanishes. On the other hand in the continuous spectrum re- 
gion G is finite, nonvanishing and complex on the real axis: as a consequence 
the same occurs for E. 

In conclusion E, and consequently m*, v, t and V, are well-behaved func- 
tions in the neighborhoods of any point of the spectrum of G. However they 
are singular elsewhere on the real axis (actually between the poles of G), in 
correspondence to the discrete part of the spectrum of the nucleus, whose 
occurrence is allowed by its finite size. 
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Now let us shortly comment the dependence of the eigenvalues upon k 
and p, which is fixed by eq. (|25|) (at least for the advanced part of the 
spectrum; the generalization however is obvious): 



ei{k, p) = ^p2 + Ml - V(P - k)2 + (M^i)2 . (121) 

This relation suggests the choice of k"^ and {p — kY as independent variables 
for expressing m*, v, t and V, because it entails 

(p - kf = {m^^-'Y . (122) 



Now the eigenvalues satisfy eq. (|99D , which, setting f = k^ — ml, can be 
recast as follows: 

fiip-ky,k') = (123) 
which in turn can be inverted to yield 

(p - kf = ^lik') . (124) 
Combining then ( 122 ) with ( |124|) one finally gets 

(^^(P) = VP (125) 

which constrains the functional dependence of the self-energy. 

As a simple example fulfilling eq. ( |124D we consider the case v = 1 and 
if: = (which implies the absence of an effective momentum) and requiring 
in addition to have m* depending upon {p — k)"^ only. Then the eigenvalue 
equation (^) reduces to 

{p - k f = m*\{p - kY) . (126) 

In this simple example the only unconstrained quantity left out is V, which 
can then be fixed by comparing with the experimental data. 

In conclusion the purpose of this subsection has been to point out how 
the structure of the self-energy is constrained, in a Lorentz covariant trasla- 
tionally invariant framework for a finite system, by the analytic properties it 
must display and by its dependence upon the momentum p of the nucleus. 
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5 Conclusions 



In this paper we propose a relativistic, hamiltonian theory of nuclei, which 
also fulfills translational invariance. 

In order to set up such a framework we focused our attention on the 
spectral amplitude: indeed by one side the latter represents a most suitable 
tool for extending a basic concept of the nuclear physics in the non-relativistic 
domain, namely the one of single particle orbital, to the relativistic regime, 
on the other side it turns out to be natural to express the energy of the 
nucleus and its 3- and 4-momentum distributions, in a scheme respecting 
covariance and in which the number of nucleons is not conserved, precisely 
in terms of the spectral amplitude. 

We achieve our scope within a model, namely assuming the nucleus to 
be a composite system of nucleons and mesons. Although the latter are not 
elementary degrees of freedom, yet we believe that they can be considered 
as "effective" building blocks useful in describing the nuclear structure in a 
quite wide kinematical region: therefore to them we have first applied our 
formalism. On the other hand if and when more fundamental degrees of 
freedom, namely the quarks and gluons, should be accounted for, then our 
theory can equally well be applied to this new dynamical situation. 

As an aside, in sect. 3 of the present paper we have also discussed the 
limiting case, opposite to the above one, where the energies involved are so 
small that the nucleus lends itself to be described solely in terms of (suitably 
renormalized) nucleonic degrees of freedom. This is the extreme situation 
corresponding to the more traditional nuclear physics where the mesonic 
degrees of freedom are eliminated in favour of static forces acting between 
the nucleons. Of course the limits of such a scheme are met when the Meson 
Exchange Currents (MEC) intervene in a decisive manner, as, e.g., in the 
electro disintegration of the deuteron at large momentum transfer . 

In the present paper we have not addressed the question related to the 
global gauge invariance of the theory, a symmetry required to consistently 
account for both the currents and the forces carried by the mesons. We have 
rather centered our attention on the dynamics of the many-body system ex- 
ploring in particular the fermionic self-energy, the basic quantity entering 
into the definition of the fermion propagator. More specifically we started 
from its relativistic, translational invariant definition and by combining its 
Lehmann representation with the relativistic Dirac's equation we succeeded 
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in deriving an equation obeyed by the spectral amplitude 0; itself. This one 
can be solved once the self-energy is known and the associated eigenvalues 
directly provide the energy spectrum of the nucleus. We have actually done 
that assuming a general Poincare invariant self-energy: thus we obtained 
the spectral amplitude and the spectrum of the nucleus. Notably the lat- 
ter turned out to be a shifted Dirac's spectrum with an effective mass and 
momentum. 

In attempting a more microscopic approach, in particular to account 
for the final state interaction, we considered a Bethe-Salpeter-type equation 
whose solution happens to be closely connected with the spectral amplitude 
0/. Worth noticing is that our Bethe-Salpeter approach, for a particular 
choice of the nucleon-nucleus interaction and for infinite nuclear matter, re- 
duces to QHD. 

Finally we have explored the analytical properties of the self-energy in 
the energy variable. We have found that, in correspondence of the zeros 
of the fermion propagator, the self-energy displays simple poles, which are 
clearly connected with the bound states of the discrete part of the nuclear 
spectrum. The locations of the poles is dependent upon the momentum 
p of the nucleus, which is perfectly defined in our translational invariant 
framework. The possibility is thus offered of exploring how the spectrum of 
the nucleus modifies as a function of p: this seems to us a fascinating topic to 
investigate, in particular we conjecture that in the infinite momentum frame 
the nuclear spectrum will considerably simplify. 

In conclusion in the present paper we have outlined a theoretical frame- 
work suitable for handling relativistic finite nuclear system. The next step 
would be to deal with concrete applications: for example we can use our 
scheme to ascertain the antinucleonic components in the nuclear wave func- 
tion, which has never been properly assessed. Furthermore, as already men- 
tioned, there is no reason preventing the extension of the present formalism 
to describe a system made by quarks and gluons, because the constrains we 
obtain on the self-energy only stem from Lorentz and translational invari- 
ance. In fact the structure of a quark self-energy in a nucleon should be 
simpler, because the excitation spectrum of the nucleon is simpler than the 
one of a nucleus. 

Clearly the deep inelastic scattering on both the nucleon and the nucleus 
offer the best testing ground of the present formalism, which might provide 
a help in achieving a better understanding of those processes. 
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A Appendix: Boosting the spectral ampli- 
tude 

We remind the definition 

0Kk,p)-<(^-l)i,k-p|^(O)|Ap> ■ (127) 

Let A be a pure boost which transforms the state | A, p > to a state at rest, 
i.e. 

A|Ap>=|Ap = 0> . (128) 

Obviously 

P) =< {A - k - p|A-^ (a^(O)A-^) A| a P > . (129) 
We also know that 

AV'(0)A-^ = ^^^(0) . (130) 

where 

S = e^-^/^)'^"'^^"'' (131) 

and 

C^alB^ ^ba,ll3\ ■ (132) 

the parameters of a pure boost having the form 

The matrix 5" transforming a particle at rest (with mass m) into one 
moving with momentum q is 



lE + m ( 1 ^ 

E+m 

with E — \/q^ + m^. The velocity of the moving particle is 

V = ^ = — 

In our case the boost is characterized by the velocity 



(134) 



P 



(135) 
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It is now immediate to conclude that 

Am^-' = Si^P)^ (136) 

Finally we need A\{A — l)i,p — 'k>. It clearly transforms in another state 
having the same internal quantum numbers but a transformed momentum, 
namely 



We define in term of q the quantity 

0(q)=<(A-l)^,q|V^(O)|AO> (138) 

and the previously derived transformation properties provide us the factor- 
ization 

j^p) 0(q) (139) 
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